clear;clc;close all
alpha_f=deg2rad(50);%напревление максимума антенны ДН антенны
beta_f=deg2rad(0);


alpha_s=deg2rad(-90:5:90);
beta_s=deg2rad(-90:5:90);%otkyda prihodit signal

beta = H (alpha_f,beta_f);%kyda napravlena antenna.
for a=1:length(alpha_s)
    for b=1:length(beta_s)
        S = 1;
        y = H(alpha_s(a),beta_s(b))*S;% вектор фокусировки
        U = beta'* y; %vihod antennoi reshetki
        F(a,b) = abs(U/S)^2; %DN
    end;
end;


[alpha_sm,beta_sm] = meshgrid(alpha_s,beta_s);%формируем поверхность

[Fx,Fy,Fz] = sph2cart(beta_sm, alpha_sm, F);


figure;
surf(Fx,Fy,Fz)
